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ABSTRACT 

We have fitted the surface brightness of a sample of 79 elliptical galaxies pertaining 
to the Coma cluster of galaxies using the Sersic profile. This model is defined through 
three primary parameters: scale length (a), intensity (So)i and a shape parameter 
{v)] physical and astrophysical quantities may be computed from these parameters. 
We show that correlations are stronger among primary parameters than the classical 
astrophysical ones. In particular, the galaxies follow a high correlation in v and a 
parameters. We show that the and a correlation satisfies a constant specific entropy 
condition. We propose to use this entropy relation as distance indicator for clusters. 

Key words: galaxies: clusters: Coma - galaxies: distances - galaxies: fundamental 
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1 INTRODUCTION 

The photometrical properties of elliptical galaxies have al¬ 
ready been extensively analysed, with special attention paid 
to the entire range of sizes and luminosities, from giant to 
dwarf ellipticals. Profile laws are in general defined by vari¬ 
ous parameters, for example, scale length (a), intensity (Eq), 
and one or more parameters specifying the shape (structure 
parameters). These parameters can be categorized as basic 
or primary, in the sense that they come out naturally from 
a mathematical definition. 

The de Vaucouleurs profile has no structure parameter; 
for a long time it appeared to be in very good agreement with 
observations, although hints of systematic deviations were 
observed; those deviations appear similar for galaxies with 
the same luminosities (Michard 1985, Schombert J.M. 1986). 
New observations, especiall y high resolution data from the 
HST (e.g. Ferrarese et al. (1994)), have definitively shown 


surface brigh tness laws. Among many possibilities, Sersic’s 
recipe (1968), 


that the de Vaucouleurs’ profile - rather too inflexible - does 
not allow a convenient description of the actual structure of 
all elliptical galaxies. 

Some attempts have been made to obtain more general 


* Based on observations collected at the Canada-France-Hawaii 
telescope, operated by the National Research Council of Canada, 
the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii 


E(r) = So exp(—u‘') , u = r/a , 


( 1 ) 


a generalization of the de Vaucouleurs-law, seems to be bet- 
t er su ited to describ e ellip tical brightness pr ofiles (e.g. Ciotti 
(1991 ) , Cao n et al. (1993), Graham et al. (1996), Courteau 
et al. ( 1996 )). 

In this letter we show the results of the photometri¬ 
cal analysis of a sample of elliptical galaxies in the Coma 
cluster using a Sersic-model (or a v model). From the basic 
parameters [a. So, u], it is possible to build physical and as¬ 
trophysical quantities (effective radius, effective luminosity, 
total magnitude, energy, entropy, etc...). All these quan¬ 
tities are in fact a combination of primary parameters; we 
therefore propose to use primary parameters instead of the 
classical ‘astrophysical’ quantities in order to obtain better 
defined correlations. 

Any correlation between a dimensionless parameter 
and a distance-dependent one has traditionally led to 
the definition of a distance indicator; for instance, the 
Cepheid luminosity-period relation (Luminosity/Period), 
the Tully-Fisher relation (Luminosity/Rotation), and the 
Faber-Jackson relation (Luminosity/Vel ocity d ispersion) are 
widely used. Recently, Young & Currie (1994) have pointed 
out a correlation between the (distance-dependent) intrinsic 
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luminosity and the (dimensionless) shape parameter of the 
Sersic-law for dwarf ellipticals. 

A theoretical understanding of the physical underlying 
processes is crucial to assure the validity of a correlation 
between the parameters describing an astronomical object. 
Here, we show that elliptical galaxies are located in the lo¬ 
cus dehned by constant specihc entropy s, in the primary 
parameters plane [v, o]. The fact that the elliptical galaxies 
in Comamainly have the same specihc entropy should some¬ 
what reflect the intrinsic properties of self-gravitating sys¬ 
tems formed through violent relaxation processes. We will 
discuss whether the entropy relation connecting v and a may 
be used as a distance indicator. 


2 THE SERSIC-MODEL 


L*{v) = 27rir(^) 


(4) 


where r(a:) is the Gamma function. L*{v) depends only on 
the structure parameter u. The total magnitude is mtot = 
—2.51ogl/tot+ cte. We have approximated the magnitude 
m*{y) = —2.5 log by the following expression: 


rn{v) = 2.199 - 1.246 v - 


2.720 

V 


0.2195 


(5) 


which is accurate to 3% for 0.1 < v < 1.0. 

• for the effective Radius: = aUe, where U^{v) de¬ 

pends only on v, we have approximated it by an analytic 
expression: 


\nUe{v) = i(0.5434 - 1.069 Ini/) (6) 

• for the effective brightness, /Te: 


2.1 Scale and structure parameters 


Surface brightness functions generally used to describe el¬ 
liptical galaxies can be written as: 


E = E.E(|d|;(u)),u = (r/a) 


( 2 ) 


were F does not depend on the normalization factor, Es. Es 
is a characteristic intensity and a is a characteristic scale. 
Both are scaling parameters, while the set a, contains the 
structural parameters which account for the specific shape 
of the surface brightness prohle. The so-called / 3-model, 7 - 
model, or more generally the {a, (3, yj-model (Zhao 1996), 
are known examples. 

From both scaling and structural parameters which are 
primary in the sense of being directly obtained from the 
ht, total luminosity, effective radius, effective magnitude, 
etc... may be calculated. Since they result from a combi¬ 
nation of the primary parameters, they can be considered 
as secondary quantities. 


/ie = -2.51og(Eo)-bl.085717e(!z) (7) 

As it can be seen, fj.^ is the sum of two terms; the first one 
is a function of Eq and the second one only depends on v. 

• For the mean effective brightness, < p >e' 

< p >e= 2.5 log 27r — 2.5 log Eq -|- 5 log Ue{v) -\- m*{u){8) 

< p >e is the sum of 4 terms, the first being a function of 
Eo while the second and third depend only on v. 

Therefore, the photometrical properties of elliptical 
galaxies can be described in terms of a combination of pri¬ 
mary parameters. That is to say, we may transform the 
triplet [Eo,a, iz] into the triplet [Ltot{mtot), Rcv]- Correla¬ 
tions between the secondary parameters [I/tot(mtot), Re, i^] 
should therefore show larger dispersions than those corre¬ 
sponding to [Eo,a,zz]. Notice that any strong correlation 
between either Eq or a, and n should provide a distance 
indicator. 


2.2 The Sersic distribution 

To proceed with the ht we have decided to use a general 
brightness prohle pertaining to the family described by Eq. 
(^), namely the Sersic model. The htting function, that we 
will refer as the i/-model (due to the analogy with the (3- 
model) is given by Eq. (^. 

The relation between v and n (which is often used) is 
n = Ijv. For v = 0.25 (i.e., n = 4) the de Vaucouleurs- 
model is recovered. We note that the de Vaucouleurs-model 


3 ANALYSIS OF THE DATA AND 
CORRELATIONS 

3.1 The data 

The sample of elliptical galaxies has been selected from the 


catalogue by Lobo et al. (1997), based on V-band CCD im¬ 


ages covering a region of 30' x 21'.6 centered on the position 
a = 12'‘59"*42'!71, <5 = -b27°58' 14'12, observed at the 3.6 m 
Canada-France-Hawaii Telescope in 1993. We have selected 
79 galaxies, which cover the range [13-17.8] in V-magnitudes 


has no structural parameter, implying that any galaxy pro- 

and [0.6-20 kpc] in ehective radius. We assume a 

distance 

hie can be deduced from a unique galaxy model by scaling 

for Coma of 137 hG Mpc ( 

Colless & Dunn 1996 

). A full 

(homology property). 

description of the sample will be given elsewhere (t 

lerbal et 


al. 1997). 


2.3 Total and effective quantities 

To calculate the total luminosity and ‘ehective’ quantities 
(i.e. quantities related to the ‘effective’ radius Re) we have 
used the classical dehnitions and, after straightforward cal¬ 
culations, we find: 

• for the total luminosity, Ltot: 

Ltot = Eoa^L*(iz) (3) 

with: 


3.2 The fitting method 

We have htted the integrated flux (the luminosity growth 
curve) as a function of the surface. We decided to use such 
a fit since the irregularities of the galaxy are smoothed out. 
While htting the light prohle is very sensitive to the outer 
parts, htting the integrated hux gives a better represe ntati on 
of the global light distribution (Prugniel & Simien (1996)). 
Moreover, in this way we can take into account the ellipticity 
of the galaxies. 
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Figure 1. Distribution of the primary parameters a, Sq (top), 
and u, So (bottom), is a dimensionless parameter, a is expressed 
in kpc h'^Q, and Sq in erg arcsec“^ 



V 


Figure 2. Distribution of the primary parameters a, ly. The units 
are as in Fig. 1. The solid line represents a constant specific en¬ 
tropy track (cf. section 4.1). 


Table 1. Correlation between the parameters of the Sersic-law. 
Both p and r are defined in the interval [—1,1], the 0 value mean¬ 
ing no correlation. Higher absolute values of 2 : indicate a greater 
significance of the correlations 


Parameters 

p-Spearman 

2 :-value 

r-Kendall 

2 :-value 

u—a 

0.97 

8.53 

0.85 

11.06 

v-T,o 

-0.80 

-7.08 

-0.62 

-8.07 

a-'So 

-0.85 

-7.52 

-0.68 

-8.88 

y-Re 

-0.54 

-4.73 

-0.38 

-4.91 

y-mtot 

0.48 

4.23 

0.34 

4.42 


We therefore write equation ([^ in an integrated form: 

2b (2 


= = j.(|) 


( 9 ) 


where 7 ( 0 , x) is the incomplete gamma function, b is related 
to the length scale o by & = nC, and e is the surface cor¬ 
responding to the region for which the flux is higher than a 
given level E (cf. Eq.(^)). The corresponding mean radius is 
then defined as r = LtotO is the total flux within 

a given e. 

The sky level was determined in a similar way as 
Jprgensen & Franx (1994). In order to limit the influence of 
the seeing (FWHM « 0.9”) the data points are taken from 
3.0” outwards, out to a surface brightness = 24.0 (i.e. 3a 
signal above the background sky noise). The fitting process 
has been done under Interactive Data Language (IDL) sup¬ 
plying data values with their l-cr errors. These are very small 
(~ 1%, 3%, and 2% for u, a, Eq, respectively) and are cor¬ 
related with each other as in the de Vaucouleurs or /3-model 
case. A more detailed d escription of the fi tting procedure 
will be given elsewhere (Gerbal et al. 1997), together with 


the resulting values for primary and secondary parameters. 


3.3 Correlations between parameters 

In Fig. 1^ and Fig. we show the relations between the 3 
primary parameters: [i/, Eq], [o, Eq], and [a,v]. 

In order to quantify the correlations in Figs. and ^ 
we have performed the Spearman and Kendall rank correla¬ 


tion (non-parametric) tests ( Siggel fc Castellan 1989 ). The 
results are presented in Table 

It can be seen that there are correlations between all 
primary parameters. However the correlations related to the 
normalization factor Eq are less pronounced, while a and u 
are strongly correlated. The greater a, the greater 1 /, i.e., the 
greater the characteristic length, the steeper the brightness 
profile. 

In Fig. ^ we show the correlations between the sec¬ 
ondary parameters \v,Re\ and [u,mtot]- The secondary pa¬ 
rameters effective radius, Re and total magnitude, nitot, are 
clearly much less strongly correlated with v. 

This result could be explained in terms of the noisy 
combination of the primary parameters to obtain the sec¬ 
ond ones and, to a certain extent, the role played by Eq 
( see eq. ^), which is not so well correlated to the two other 
primary parameters. 


4 IMPLICATIONS OF THE RELATION [a, u] 
4.1 A physical interpretation 


As noted by Graham et al. (1996), the physical origin of the 
[a, f] correlation is not understood up to now. We will show 
that this relationship has a simple physical interpretation 
in terms of gravo-thermodynamics. Notice that this kind of 
relation exists for a large family of obje cts - from the nuclei 
of spiral galaxies (Courteau et al. 1996) to Brighter Cluster 
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Figure 3. Correlations between u and two secondary parameters: 
Re (in kpc h'^Q , top) and mtot (in magnitude, bottom). 


gala xies, and even for t he X-ray emitting gas in galaxy clus¬ 
ters (Durret et al. 1994). Therefore, the physical explanation 
should be only related to a and v and to a small extent to 
So. 


Gravitational dynamical systems in quasi-equilibrium 
are in a state of quasi-stationary entropy. We study here 
the variation of the entropy with the parameters a and u 
of the Sersic-law. We use a modified deprojection of the de 
Vaucouleurs law derived by Mellier and Mathez (1987) to 
represent the 3D-density: 


p{r)^po{^ %xp(-Q‘^) , (10) 

where po is the normalization parameter. For any Vi there is 
one Pi so that the r'i-model should be recovered by project¬ 
ing p(r) in 2-D. When v = 0.25 then p = —0.855. In order 
to recover the Sersic-law from the projection of Eq. (0) we 
have derived, following Mellier & Mathez, the relation: 

p{u) = 0.9976 - 0.5772;/ + 0.03243;/^ , (11) 


which is accurate to better that 0.1% in the range 0.1 < v < 

1 . 0 . 

Assuming spherical symmetry, hydrostatic equilibrium, 
and an isotropic velocity tensor, it is straightforward to com¬ 
pute the total mass, Mtot (assuming a constant mass-to- 
light ratio), potential energy, U, and pressure profile, P{r). 
Finally, the entropy (5), and t he specific entropy (s) (fol¬ 
lowing White & Narayan (1987), assuming an ideal gas), are 
given by: 

^ [ ln(p-®/^p3/^)pdy , (12) 

Mtot Jy 


Mu 


the integral is over all the volume V of the galaxy. 

In Fig. ^ we have plotted the specific entropy as a func¬ 
tion of a and u. As can be expected, the entropy of a self- 
gravitational system of unlimited extent has no maximum. 
With the present calculation it increases with u and with 
decreasing a. Superimposed to the specihc entropy curve 
(isentropic track) are the observed values of a and v from 
the fits; one can note that they closely follow one another. 

This result suggests that there should exist some in¬ 
trinsic relation between a and u, the structural parameters 
of the Sersic-law. However, in order to be able to obtain 
a maximum entropy in a dynamical system, one needs to 
impose more constraints, in addition to Eq. ([l^). The isen¬ 
tropic track occupied by the observed values of a and i/ may 
then be understood as a maximum entropy class of some 
more restricted distribution. 


4.2 Towards a distance indicator 

The ;/-model, contrary to the de Vaucouleurs-law, does not 
allow to deduce any brightness prohle just by scaling; the ho¬ 
mology property of the de Vaucouleurs-law is lost. Instead, 
the shape property (v) is related to the scale (a) by the re¬ 
lation of constant entropy (^|), s = so- Since a depends on 
the distance whereas u does not, their correlation can be 
used as a distance indicator. 

All the galaxies used for the ht belong to the Coma 
cluster, so they are at the same distance (the cluster mem¬ 
bership has been determined with the redshift, see Biviano 
et al. 1995), and thus provide us with a calibration for the 
distance indicator based on the a-u relation. If the entropy 
relationship were ‘universal’, i.e., if it were the same for dif¬ 
ferent clusters, it could be used as a robust distance indica¬ 
tor. To any observed v^k) of any galaxy gfe in the Universe, 
would then correspond an ak (in angle units) which may 
be compared to the corresponding value obtained from the 
entropy-relation (Fig. 2), allowing to derive the distance to 
gk- 

One could wonder if the a-u relation is weakened by 
the correlation between the fitting parameters for a given 
galaxy. Since the isentropic track is predicted by theory, we 
actually believe that the a-v relation is robust. 

The dispersion of the residue distribution of A = 
(ttcaic — n) is — 0.1 leading to a 1-cr error of 10% of the 
residue. The observed dispersion of A could be due to some 
difficulties in the htting process, or due to physical rea¬ 
sons. In fact, it seems reasonable to imagine that a galaxy 
which has undergone some energy and entropy exchanges 
with other galaxies (merging, tidal effects, etc...) could not 
verify the entropy relation. 

Notice that the z/-profile is a 2-D brightness profile, 
whereas a 3-D density profile is needed to calculate the spe¬ 
cific entropy. We have used a constant mass-to-light ratio 
to transform from one to the other. Therefore, those galax¬ 
ies being more luminous than normal ones (for instance a 
galaxy hosting a starburst), will be out of the theoretical 
relationship. 
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5 CONCLUSIONS 

We have given a physical interpretation for the observed 
relation between two of the primary parameters of the v- 
profile: the values of the a and u parameters of Coma cluster 
galaxies approximately have a constant specihc entropy. The 
dispersion around the best fit is small, therefore allowing the 
use of this relation as a distance indicator (relative to the 
Coma cluster). 

We must stress that both the observational and theoret¬ 
ical evidences found in this work have to be verified. On the 
one hand, it is necessary to quantify the universality (i.e. in 
other clusters) of the (a, u) correlation found for the galaxies 
in the Coma cluster. On the other hand, the validity of the 
assumptions made to derive the theoretical entropy relation¬ 
ship (constant M/L, ideal gas entropy) must be confirmed. 
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